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We suggest that the lattice Dirac spectra in QCD at fi- 
nite temperature may be understood using a gaussian unitary 
ensemble for Wilson fermions, and a chiral gaussian unitary 
ensemble for Kogut-Susskind fermions. For Kogut-Susskind 
fermions, the lattice results by the Columbia group are in 
good agreement with the spectral distribution following from 
a cubic equation, both for the valence quark distribution and 
the anomalous symmetry breaking. We explicitly construct 
a number of Dirac spectra for Wilson fermions at finite tem- 
perature, and use the end-point singularities to derive ana- 
lytically the pertinent critical lines. For the physical current 
masses, the matrix model shows a transition from a delocal- 
ized phase at low temperature, to a localized phase at high 
temperature. The localization is over the thermal wavelength 
of the quark modes. For heavier masses, the spectral distri- 
bution reflects on localized states with competitive effects be- 
tween the quark Compton wavelength and the thermal wave- 
length. Some further suggestions for lattice simulations are 
made. 



1. Numerous lattice simulations indicate that massive 
QCD undergoes either a rapid cross-over or a phase tran- 
sition at high temperature. The character and nature of 
these phase changes are still debated M. Most numeri- 
cal analyses on the lattice have relied on Kogut-Susskind 
(staggered) quarks, although results using standard or 
modified Wilson actions are now becoming available [Q. 

QCD with three flavors undergoes a first order tran- 
sition for light current quark masses. The transition is 
believed to switch to a second order for zero up and down 
quark masses with a heavy strange quark, and a cross- 
over for all other finite quark masses. When the quarks 
become very heavy, they decouple. A pure Yang-Mills 
phase undergoes a first order transition at high tempera- 
ture. The precise values of the critical temperature, and 
critical quark masses are not yet unambiguous, as they 
seem to depend on scaling and the way the quarks are 
set on the lattice || . 

The lattice results have spurred a number of theoreti- 
cal investigations using mostly effective models, all aimed 
at understanding the whys behind the disappearance of 
chiral symmetry whether through a sharp phase transi- 
tion or a smooth cross-over B. Most noteworthy are 



the suggestions made by Pisarski and Wilczek concern- 
ing the generic character of the chiral phase transition 
and its relation to scaling and universality ||. 

The purpose of this letter is to use some insights from 
random matrix theory || , to analyze some aspects of chi- 
ral symmetry in QCD at finite temperature and nonzero 
current quark masses [@-||] . In section 2, we show that the 
recent finite temperature results by the Columbia group 
ifTTf using two-flavor QCD with staggered fermions, in 
the semi-quenched approximation are well described by 
a high temperature version of the chiral gaussian unitary 
ensemble with thermal masses, as expected from dimen- 
sional reduction. The spectral distribution follows from 
a cubic (Cardano) equation. The issue related to the lat- 
tice simulation of the anomalous symmetry breaking is 
also recovered, although the chiral random matrix model 
does not account for the U(l) anomaly. In section 3, we 
suggest that three-flavor QCD with Wilson fermions, in 
the quenched approximation, follows from a gaussian uni- 
tary ensemble. The role of the various Matsubara modes 
is highlighted. In section 4, the Dirac spectra for different 
masses and temperature are just given by a superposition 
of three solutions to a cubic (Cardano) equation. The 
phase diagram of the random matrix model is shown to 
follow from the end-point of the Dirac spectrum analyt- 
ically. Our conclusions and recommendations are given 
in section 5. 

2. Recent lattice simulations by the Columbia 
group pd| ] using staggered fermions have unraveled inter- 
esting aspects of the Dirac spectrum of two-flavor QCD 
in the semi-quenched approximation on a 16 3 x 4 and 
32 3 x 4 lattices. If we were to denote by £ the valence 
quarks with mass m^, then the valence quark condensate 
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1 g(\) ~ |A| 3 asymptotically in QCD, and ([l]) diverges 
quadratically. The finite lattice spacing provides a natural 
ultraviolet cutoff, making ([[]) finite. A better definition of ([j]) 
is on a cooled lattice. In random matrix models, the spectra 
are bounded. 
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where g(X) is the Dirac eigenvalue distributions associ- 
ated to ip for two-flavor and massive QCD. The third 
flavor Q is not included while averaging over the sea 
fermions. The dependence on the sea mass m s (not to 
be confused with strangeness) and the number of fla- 
vors in g(X), stem from the random averaging over the 
gauge configurations in the presence of the two-flavor 
fermion determinant ]fl[ |. For two degenerate flavors 
with m s a — O.Of, the behavior of (jl]) versus on a 
f6 3 x 4 lattice is shown in Fig. 1, for eight temperatures 
(3 = 6/g 2 around the critical temperature (3 C = 5.275. 
The transition in this case is believed to be second order, 
with a critical temperature T c — l/4a ~ 150 MeV 0. In 
physical units, the lattice spacing at the critical temper- 
ature is a ~ .33 fm, and the sea quark mass is m s = 6 
MeV. 
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FIG. 1. Semi-quenched condensate versus the valence 
quark mass for two-flavor QCD Jll| with a sea mass 
m s a = 0.01. 

The bulk of the staggered lattice results can be under- 
stood from simple random matrix theory [@-^) , if we note 
that for a bounded spectral density, the valence conden- 
sate is just 
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which is just the spectral density for A = irri(. For stag- 
gered QCD fermions at zero temperature, the spectral 
distribution can be mapped on a chiral gaussian unitary 
ensemble, provided that the gauge configurations are suf- 
ficiently random. At finite temperature, this is equivalent 
to an assembly of Matsubara modes, each described by a 
gaussian unitary ensemble. At high temperature the only 
relevant modes are ±7rT, as suggested by the screening 



length measurements [|| , and dimensional reduction jl2) . 
We note that ttT is the minimum wavenumber k = 2ir/X 
of a massless quark mode on a circle of length 1/T with 
a wavelength A = 2/T. For massive quarks, the relevant 
modes are ±M(T) = ± v / m 2 + tt 2 T 2 @. With this in 
mind, a suggestive high temperature chiral random ma- 
trix model for two degenerate and massive flavors is 
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where each entry in (||) is N x N valued. The i?'s are 
described by a gaussian unitary ensemble. The off-block- 
diagonal structure of (||) is suggestive of the chiral struc- 
ture of ip, which is preserved in the staggered fermion 
formulation. This is not the case for the Wilson fermions, 
as we discuss below. All the scales in (|3|) are given in units 
of 1/S = (100 MeV)" 1 . Although (g) is well-motivated 
at high temperature, we will use it to investigate the tem- 
perature ranges around the critical temperature. This is 
a bit justified by the fact that the lattice simulations 
of the spatial screening lengths show a behavior that is 
consistent with the occurrence of the lowest Matsubara 
modes around T ~ T c . For m s = 0, and in the case of 
one-flavor and a single Matsubara mode, (^) has been 
numerically investigated in j§] . 

The spectral distribution associated to (||) follows from 
the discontinuity of the pertinent solution to a cubic 
(Cardano) equation, 

G 3 (z) - 2zG 2 {z) + {z 2 - M 2 (T) + 1)G - z = (4) 

with — Trg(X) = ImG(A + i0), this was also observed in || 
for the massless case. The equation (||) follows readily 
from the law of addition of random plus deterministic 
matrices [10[. 

In Fig. 2Twe show the behaviors of the valence quark 
condensate for two choices of the sea quark masses, m s — 
and 0.1, and various (temperature) values as obtained 
from the Cardano solution (0). In physical units, m = 
and 10 MeV. We have identified (3 with T, through^] 



n(T - T c ) = [fl - (3 C ) 



(5) 



with the critical temperature T c in the chiral random ma- 
trix model and (3 C = 5.275 as suggested by the lattice cal- 
culations |]ll"f . The transition in the chiral random matrix 
model is mean-field in character (large N). For m s = 0.1 
(10 MeV) the random matrix model seems to follow qual- 
itatively well the lattice results for m s a = 0.01 (6 MeV), 



2 Note that here (3 is explicitly related to the temperature, 
while on the lattice (3 = 6/g 2 is implicitly related to the tem- 
perature through the coupling constant after a change in the 
lattice spacing. The use of the same parameter for both de- 
nominations is only suggestive. 
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except for < 10~ 5 , where the spectrum becomes sen- 
sitive to the finite size of the lattice, as illustrated by 
the bending of the upper curves. Finite size effects set in 
when the pion Compton wavelength becomes comparable 
to the lattice size. Using the PCAC relation this implies 
that m c < CC >> l/iV c y 4 . For V 4 = 4 x 16 3 , this puts 
a lower bound on the valence mass > 10 , which is 
about consistent with the lattice results. We note that 
the above interpretation of the finite temperature lattice 
simulation is different from the zero temperature univer- 
sality argument used in jl3| ] , but in qualitative agreement 
with the finite-N numerical analysis carried in S , in the 
massless case. 

m =0 
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FIG. 2. Semi-quenched condensate versus the valence 
quark mass for two- flavor QCD, obtained from eq.(^) 
with the Cardano solution of eq.(W), for the sea quark masses 
m s — (upper) and m a =0.1 (lower). 



depends sensitively on the value of the sea quark mass, as 
can be seen by comparing the upper and lower curves on 
Fig. 0. For m s = 0, the critical temperature is (3 = 5.275 
(fifth line from the top), whereas for m s = 0.1, the crit- 
ical value is (3 — 5.27 (the fourth line from the top). In 
the random matrix model (Cardano solution), the depen- 
dence of the critical temperature on the sea quark mass is 
simply A/3 C = y/l — m 2 s . Both lines have the same crit- 
ical exponent 8 = 3 (mean- field) . This result was first 



-,0.33 



For 



established in g|. For (3 = 5.275, < CC >'~ <<\ 
(3 = 5.34, < CC >— 8™^. These values are to be con- 
trasted with the ones quoted by Chandrasekharan and 
Christ O), namely < CC >~ ™c' e f° r P = 5.275, and 
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.0.98 



for (3 — 5.34. The closeness of the latter 



to the random matrix result, suggests perhaps an un- 
derestimation of the critical exponent 5 = 5/3 (ll]] in 
comparison with the mean-field one 5 = 3, possibly due 
to the large finite sizes effects mentioned above as visi- 
ble in Fig. [l] (dropping plateau's for small masses). We 
recall that a two-flavor simulation of QCD using finite- 
temperature cumulants yields 8 in the range 4-5 [fL4| , 
hence closer to the mean- field result. 

The fact that the semi-quenched condensate appears 
to vanish linearly with the valence quark mass rriQ in the 
high temperature phase ((3 > (3 C ) prompted the Columbia 
group to analyze the asymmetry 
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between the isotriplet and its axial partner, 
bounded spectrum Q this expression reduces to 



Trg(im) d 
m dm 



(6) 
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Below the critical temperature and as m goes to zero, 
the first term diverges (pion pole), while the second term 
approaches zero (plateau), provided that g(im) does not 
develop a low fractional dependence on the sea quark 
mass (set equal to the valence one). Above the critical 
temperature, the first term vanishes for any value of the 
mass m, which is followed by an infinite jump in oj at the 
critical temperature. For finite values of the mass m this 
jump is smeared out. 

Figure ||a shows the behavior for the Cardano solu- 
tion with m = 0.005 (0.5 MeV). Figure |b shows the 
same curve for ma = 0.01 (6 MeV) as obtained by Chan- 
drasekharan and Christ on a 16 3 x 4 lattice with stag- 
gered fermions pd| ]. The similarity is striking, although 
the mass scales are different. The fact that in the ran- 
dom matrix model, g{im), does not develop a low frac- 
tional dependence on m is traced back to the structure 
of the thermal masses (T-regulated) . We recall that the 



We note that the reading of the critical temperature f3 c 



See footnote 1 above. 
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random matrix model considered here has no bearing on 
the axial-singlet anomaly. Its discussion in the context 
of random matrix models goes outside the scope of this 
work. 
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FIG. 3. (a) Results from the random matrix model for 
m = 0.005 in units of E = 1. (b) results for staggered fermions 
on a 16 3 x 4 lattice with ma = 0.01 Ol. 



3. Wilson fcrmion formulations are more natural to 
use for three-flavor QCD. Recently Kalkreuter has per- 
formed detailed lattice simulations with Wilson fermions 
for two-color and two-flavor QCD |fT5| . His results are in 
qualitative agreement with random matrix theory [ ^6|Jl0| 
for the bulk spectrum. The issue with Wilson fermions, 
however, is that the Wilson terms needed to eliminate 
the lattice doublers break explicitly chiral symmetry. It is 
however expected, that in the continuum limit the break- 
ing is mild enough (irrelevant operators) to allow for a 
full restoration of chiral symmetry, modulo mass terms. 

In terms of chiral random matrix models, the perti- 
nent ensemble for investigating Wilson fermions for QCD 
is the gaussian unitary ensemble J], provided that the 
quarks are in the fundamental representation. In the 
presence of Wilson r-terms, the Dirac operator is only 
hermitean. It is not block-off-diagonal (chiral) . Random- 
ness, however, is sufficient to cause the spectral distribu- 
tion to develop a non-vanishing value at zero virtuality 
and in the chiral limit. So, although the spectrum is not 
paired, the distribution is even. This is easily understood 
if the diagonal disorder (r-terms) is not large enough [jl0| . 
In fact, this can be easily investigated by considering a 
linear combination of two chiral random matrix models. 
One with off-diagonal randomness (chiral) plus one with 
diagonal randomness (r-terms) Jl8|. 

For three flavors, we define the hermitean Dirac oper- 
ator Q = 0/75 TO/), with / = u,d,s. Here Q is 
the Dirac Hamiltonian in five-dimensions. The Dirac 



4 For QCD with two colors the relevant ensemble is the gaus- 
sian orthogonal ensemble JlTj ], 



matrix 75 plays the role of /?-Dirac matrix in (3+1) di- 
mensions. On a cylinder, the pertinent random matrix 
model for spatially constant Wilson modes is 
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where d T is the r-derivative on a cylinder of length 
P = 1/T, and R a r-independent random matrix with a 
block structure specified by eq.(|i~0|) below. Equation (jg) 
is defined with anti-periodic boundary conditions. It is 
diagonal in frequency space using 



V>(t) 



(9) 



where u) n = (2n+l)irT are the usual Matsubara frequen- 
cies. Here i{j(t) is an iV-dimensional vector for each of 
the three flavors: up, down and strange. For each Mat- 
subara frequency, and for each flavor (g) describes a set 
of N fermions, with a hermitean Hamiltonian 



Q 



lUJ n 
-TO/ 



R 



f 



(10) 



This is a the sum of a deterministic and a random piece. 
Each entry in the first matrix is N x N valued. The 
random matrix is 2N x 2N valued. In terms of (|ic|), 
eq. (0) is block-diagonal in frequency and in flavor space, 
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The random distributions for each Matsubara mode 
are decoupled, each with a weight 



P(Rr, 



1 



exp[-2NTrV(R n )} 



(12) 



This is a direct consequence of our assumption that R in 
(||) is r-independent. If we were to relax this condition, 
then the random matrix may cause the Matsubara modes 
to couple jig] . In this paper, V will be chosen gaussian, 
although polynomial weights are also possible pTgj ] . The 
resolvent of ( |ll| ) is 
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The averaging is carried with the distribution (|T^). In 
terms of ( |l3|) the spectral distribution associated to (|sj> 
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is related to the discontinuity of G{z) through the real 



i/(A) = --ImG(z = X + iO) 
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(15) 
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Equation ([li]) involves the sum over a large number of 
Matsubara modes, with N* — J2 n - 

Using either the law of addition of random matrices 
ful^of or diagrammatic techniques j2M, it follows that 
in large N the resolvent G(z) is 



G(z) = tAtY G{(z) 
v ; NfN* ^ ™ v ' 

J n,f 



(16) 



where each of the G^(z) satisfying Pastur's equation |25 



G f n (z) = ^ 
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z-Gliz)- Mi z-GL{z) + Mi 



(17) 



is a 2 x 2 diagonal matrix, with M£ (T) = 

At high temperature, only few modes contribute to 
(|l6|), i.e. luq = u—i = 7rT. Taking m u ~ rrid, the problem 
reduces to 



G(z) = -G™(z) + 



(1-a) 



G s (z) 



(18) 



with a = 1 - 1/Nf, N f = 3 and G f (z) defined by ((T^). 
( |l8| ) is a combination of three third order algebraic equa- 
tion (Cardano class) for the resolvent G(z). 

The spectral density follows from the discontinuity 
of G(z) along the imaginary axis (|l5|). In the chiral 
limit, the discontinuity is related to the chiral conden- 
sate through the Banks-Casher argument, 



< qq >= Im G(z = iO) 



(19) 



Although @ was derived in the high temperature limit, 
we expect it to be rather accurate near the critical points 
because of universality. Indeed, at the critical point and 
in large N only symmetry is relevant p3f . Below, we will 
provide an analysis of the spectral distribution associated 
to @. 

4. The explicit solution of the third order equation 
is available in algebraic form, and was discussed in flio| |. 
The combination of several third order equations allows 
for a rich phase structure. In the case of two light and 
one heavy flavors, the spectrum should exhibit at most 
four phases Pi, P2,Ps and P4. The support of the spec- 
tral function may consists of one, two, three or four dis- 
connected arcs, depending on the choice of the external 
parameters m u , m s , and T. 

Figure |^ shows the behavior of the spectral function 
for the physical choice of m u = rrid — 0.075, m s = 1.5. 
All dimensionfull quantities are measured in units of 
E = 1. In physical units, m u = nid = 7.5MeV and 
rn s = 150MeV. At zero temperature, the system is in the 
Pi phase, and the value at zero virtuality (A = 0) corre- 
sponds to a non-zero condensate. The mass of the strange 
quark is 20 times higher, therefore the spectrum feels it 
by developing two symmetric shoulders. The mass is not 



high enough, however, to provide the decoupling of the 
strange quark (P 3 phase), with symmetric and discon- 
nected humps. The effect of the temperature is to quench 
the maximum at the origin. At the critical temperature 
T*, the quarks localize, resulting into a transition from 
Pi — ► P2. For increasing T, the thermal quark modes 
contribute democratically to the spectral distribution, 
y/ml + tt 2 T 2 -> ttT, M s = y/m 2 s + tt 2 T 2 -> ttT. 



M, 




FIG. 4. Spectral function for Wilson fermions with 
m u — md = 0.075 and m s = 1.5. 




FIG. 5. Spectral function for Wilson fermions with 
m u — rrid = 0.075 and m s = 5. 

Figure || shows the same spectral distribution as Fig. |J, 
but this time for a large strange quark mass m s — 5, or 
m s = 500 MeV in physical units. The two additional 
humps at T = are just the delocalized strange quarks. 
An analysis of Cardano's solutions show that the criti- 
cal value for which this happens is 2m 2 = 11 + VT25 for 
m u = 0. Here, the initial spectrum is in the P3 phase. 
The effects of the temperature is to cause localization of 
the light quarks, first transition, from P3 to P4, resulting 
in the vanishing of the density of states at zero virtuality. 
This is followed by a second transition at higher temper- 
ature T** from P4 to P2 in which all the three flavors 
become thermally heavy. 

Figure ^ shows the same spectral distribution for three 
large but comparable masses, m u = rrid = 2 and m s = 4. 
There is no phase transition, the system is always in the 
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heavy (localized) phase P^. There is no accumulation of 
quark states at zero virtuality. The thermal effects cause 
the localization to shift from the Compton wavelength to 
the thermal wavelength. 




FIG. 6. Spectral function for Wilson fermions with 
m u — md = 2 and m s — 4. 

Figure Q shows yet another choice of heavy quark 
masses m u = = 2 and m s = 6, for which the spec- 
tral distribution undergoes a phase transition from P4 to 
P 2 . At low temperature the quarks are massive but dis- 
tinct, at high temperature the distinction is washed out 
by temperature. 




FIG. 7. Spectral function for Wilson fermions with 
m u — md ~ 2 and m s — 6. 

The spectra discussed above for the random matrix 
model, allow for an assessment of the critical lines (phase 
diagram) in terms of the end-points of the spectral distri- 
bution. The equations for the end-points are given by a 
polynomial equation that is always one degree lower than 
the one following from Pastur's equation. This is best 
described using Blue's functions pofl , i.e. functional in- 
verses of Green's functions. Generically, the Blue's func- 
tion associated to the Green's function G(z) is given by 



B{G{z)) 



(20) 



The substitution z — > B(z) and the use of ( p0|) in Pastur's 
equation (|l^) yield the pertinent equation for the Blue's 



functions. Since the Green's functions behave as \J A — z 
when z approaches the endpoint A, then 

dG _ dB_ 

~^ \end— -point — — \end— point — tJ l^4J 

For our case, this yields a quadratic equation [fill. The 
end points are explicitly given by 
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for all masses and temperatures. In the case of Figure || 
and 0, the transition from P 4 to P 2 is denned as = 
A~ . In the high mass (temperature) limit this condition 
is equivalent to M S (T*«) - M„(TU) = 2y/2. At zero 
virtuality the spectral distribution vanishes for ttT* = 

V 1 - m l- 

5. We have shown that the chiral gaussian unitary 
ensemble for two massive flavors and finite temperature, 
allows for a quantitative description of the recent simu- 
lations by Chandrasekharan and Christ using staggered 
fermions. The agreement with the random matrix model 
at p = 5.34 and higher suggests that the fit to the lattice 
data with 8 = 5/3 may be underestimating the critical 
exponent at [3 — 5.275 in comparison with the mean-field 
result by about a factor of two. We have also provided a 
qualitative understanding of the lattice anomalous sym- 
metry breaking effect in a chiral random matrix model 
without the axial-singlet anomaly, albeit with a small sea 
quark mass. 

We have suggested that in the chiral limit the gaus- 
sian unitary ensemble is the pertinent random ensemble 
for investigating Wilson fermions in QCD. We have ex- 
plicitly shown that the problem of three flavor QCD in 
the Wilson formulation is amenable to a random plus 
deterministic problem for few Matsubara modes. We 
have used arguments at high temperature to unravel a 
number of properties of the spectral distribution for the 
hermitean Dirac operator (Hamiltonian in five dimen- 
sions) . QCD with two light and one heavy flavor exhibits 
a spontaneously broken phase at low temperature, and a 
localized phase at high temperature. The localization is 
determined by the thermal wavelength. For QCD with 
massive flavors, the localization is a competitive effect 
between the thermal wavelength and the Compton wave- 
length. 

We have explicitly used the end-points of the Dirac 
spectrum to construct the critical lines for massive QCD 
at finite temperature. The critical lines (phase diagrams) 
were shown to follow from a quadratic equation fol- 
lowing from the Blue's function. Given the qualitative 
agreement noted between our results, and the staggered 
fermion calculations, we look forward to a detailed com- 
parison with future simulations using Wilson fermions 
for full QCD. We note, however, that our analysis does 
not account for the axial-singlet anomaly. Its role in the 
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context of random matrix models will be discussed else- 
where. 

The phase transitions discussed here are all mean-field 
in character. This is unavoidable in the large iV anal- 
ysis performed here. The quantitative similarities be- 
tween the Columbia results at finite temperature and 
the random matrix results, suggest that the phase tran- 
sition in quarkish QCD may be characterized by a nar- 
row Ginzburg window (perhaps of order 1/N), thereby 
asserting the lore of mean-field exponents. Within the 
Ginzburg window scaling arguments and universality 
should hold §. 

The spectral distributions we have discussed here are 
of course amenable to finite N corrections. This could be 
discussed either in the context of the corrections to Pas- 
tur's equations or using diagrammatic techniques. Also, 
specific level correlations can be discussed either in the 
staggered formulation, or the Wilson one, both for mas- 
sive QCD and finite temperature. These and other issues 
will be discussed elsewhere. 
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